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INTRODUCTION 


''The  Integrated  Refractive  Effects  Prediction  System  (IREPS)  should 
consider  sea  clutter  effects  winder  atmospheric  ducting  conditions.  The  IREPS, 
windergoing  research  and  development  at  Naval  Ocean  Systems  Center  (NOSC),  is 
intended  to  provide  a  shipboard  environmental  data  processing  and  display  cap¬ 
ability  to  assess  refractive  effects  of  the  lower  atmosphere  for  naval  EM 
systems.  Although  intended  to  be  incorporated  eventually  as  a  part  of  the 
Tactical  Environment  Support  System  (TES3),  the  IREPS  is  currently  configured 
as  an  interim  version  based  on  a  Hewlett-Packard  9845  desktop  calculator.  A 
comprehensive  discussion  of  the  IREPS  capabilities  is  presented  by  Hitney  et 
al.  (1981),  while  a  discussion  of  the  IREPS  propagation  models  is  presented  by 
Hattan  (1982). 

Numerous  limitations  of  the  IREPS  models  are  listed  by  Hitney  et  al. 
(1981),  among  which  is  the  fact  that  surface  clutter  (land  or  sea)  is  not 
included  in  the  calculation  of  radar  detection  ranges.  A  clutter  computer 
model  was  developed  and  was  first  described  by  Snyder  (1979).  Some 
modifications  have  subsequently  been  made  to  this  model  and  the  updated  model 
has  been  implemented  in  a  computer  program  in  FORTRAN  (the  language  of  the 
original  IREPS).  This  updated  model,  along  with  some  details  of  the  computer 
implementation,  is  discussed  in  this  document.  Also  presented  is  an  analysis 
indicating  the  sensitivity  of  the  model  to  variations  in  the  environmental 
parameters  required  as  input  to  the  computer  program. 

Snyder  (1979)  concluded  that  the  sea  clutter  model  developed  appeared  to 
be  adequate,  at  least  qualitatively,  for  an  elevated  layer,  ground  based  duct 
environment,  but  that  some  sort  of  quantitative  "validation"  effort  should  be 
conducted  before  implementation  into  IREPS  for  shipboard  use.  Because  of 
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manpower  and  money  limitations,  an  extensive  validation  effort  could  not  be 
undertaken.  However,  using  essentially  "free"  facilities  on  a  "not-to- inter¬ 
fere"  basis,  validation  data  were  acquired  for  a  single  location  and  a  single 
frequency.  That  validation  effort  is  also  discussed  in  this  document. 


2 


THE  PROPAGATION  AND  CLUTTER  CROSS  SECTION  MODELS 


Snyder  (1979)  and  Hitney  et  al.  (1981)  both  point  out  the  occurrence  of 
two  types  of  atmospheric  ducts  which  can  be  expected  to  greatly  influence  sea 
clutter.  One  is  the  evaporation  duct,  which  is  formed  by  the  sea  surface  and 
the  near  surface  minimum  in  the  refractive  index  prolile,  which  is  in  turn 
determined  mainly  by  the  vertical  distribution  of  water  vapor  resulting  from 
evaporation.  The  other  is  a  surface  based  duct  formed  by  an  elevated  refrac¬ 
tive  layer. 

The  propagation  environments  resulting  from  the  two  types  of  ducts  are 
uniquely  different  and  can  be  expected  to  result  in  two  different  effects  on 
sea  clutter.  The  evaporation  duct  would  tend  to  enhance  clutter  signals  over 
extended  areas,  whereas  the  ground  based  duct  formed  by  an  elevated  layer 
would  result  in  clutter  signal  enhancement  only  from  discrete  ranges. 
Propagation  in  the  two  different  environments  would  also  best  be  considered 
using  two  different  formulations.  Vertical  dimensions  for  the  evaporation 
duct  are  on  the  order  of  metres  to  very  few  tens  of  metres  and  a  waveguide 
mode  formulation  would  be  appropriate  for  this  type  of  duct.  Vertical 
dimensions  for  the  ground  based  duct  formed  by  an  elevated  refractive  layer 
are  typically  on  the  order  of  many  tens  of  metres  to  a  few  hundred  metres  and 
ray-optical  formulations  would  be  appropriate  here. 

The  clutter  model  developed  for  the  IREPS  is  based  on  a  ray-optical 
propagation  formulation  and  is  therefore  most  applicable  to  the  surface  based 
duct  formed  by  an  elevated  refractive  layer.  The  geometrical  basis  for  the 
propagation  formulation  is  developed  in  the  following  way.  The  propagation 
environment  is  assumed  to  be  two  dimensional,  with  the  earth  represented  by  a 
cylinder  and  propagation  in  a  plane  transverse  to  the  cylinder  axis.  Thus, 
earth  curvature  is  considered  in  the  direction  of  propagation  but  is  ignored 


3 


transverse  to  the  propagation  direction.  A  classical  conformal  transformation 
to  Maxwell's  equations  in.  cylindrical  form  is  employed,  converting  the  geome¬ 
try  from  cylindrical  to  rectangular  (Cartesian) .  This  transformation  has  been 
discussed  in  the  context  of  earth  flattening  by  Richter  (1966)  and  has  been 
widely  applied  in  radiowave  and  radar  propagation  problems  when  earth  curva¬ 
ture  in  the  direction  of  propagation  is  an  important  factor.  The  net  effect 
of  the  transformation  is  to  change  from  the  cylindrical  geometry  with  refrac¬ 
tive  index  given  by  n  to  the  rectangular  geometry  with  refractive  index  given 
by  n  multiplied  by  exp(z/a),  where  a  is  the  radius  of  the  earth  and,  to  first 
order  in  1/a,  and  z  is  essentially  the  altitude  above  the  aartii  (see  Figure 
1). 

The  ray  path  trajectories  are  determined  from  the  classical  eikonal 
expression  (Jones,  1979) 


,  dr 
a 

~  (n  — )  =  Vn 
ds  ds  ~ 


(1) 


where 


*• 


s  «  ^rc  length  along  a  ray  path 
n  =  n(r)  =  index  of  refraction  at  r 

A/  A/ 

AAA 

r  =  position  vector  (xr  +vr  +zr) 
~  x  J  y  z 

AAA 

x,  y,  z  =  unit  sectors 

Let 


t. 
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(i  =  x,  y,  z) 


(2) 


and  use  the  fact  that  (Vn),  =  dn/dr, 

~  i  i 


to  rewrite  ( 1 )  in  rectangular  form  as 
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Rewrite  (3a)  as 


dn  ^x  _  dn 
fcx  ds  +  n  ds  ~  dx 


(4) 


With  dn/ds  in  component  form 


dn^andx.andy.Qn^^  y  ^n_ 

ds  ~  3x  da  5yds  5z  ds  L  5r.  i 
J  x 


equation  (4)  becomes 


dtx  _  dn 
ds  5x 


(5) 


which  can  be  generalized  to  all  three  coordinate  directions.  Assuming  homo- 

a  n  b  n 

geneity  in  both  horizontal  directions  (—  =  —  =  0)  the  components  of  (5) 

dx  5y 


I'jcome 


dtx 

HI- 


t  t  j 

x  z  dn 
n  dz 


(6a) 


z  dn 
n  dz 


(6b) 
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dt 

2 

ds 


1*1  (1-t2) 
n  dz  u  rzJ 


(6c) 


If  t  =  0  initially,  corresponding  to  a  ray  having  no  initial  transverse 

component,  then  (6b)  shows  that  t  =  0  always.  Making  this  assumption  and 

transforming  to  the  angular  variable  given  by  t  =  cos  6  and  t  =  sin0  (0  is 

x  z 

the  ray  angle  with  respect  to  the  horizontal)  leads  to  the  expressions 


d0  _  1  dn  *  dz 

_  _  _  _  cos0  .  _ 


=  sin  8 


(7) 


The  refractive  index  in  (7)  is  given  in  terms  of  the  original  (cylindri¬ 
cal)  refractive  index  as 


feV  ’1 


Vs*  W* 

►  *  .  *. 


n  =  n(z)  exp(z/a) 


(8) 


Because  the  atmospheric  refractive  index  n(z)  is  only  slightly  greater  than 
unity,  it  is  convenient  to  write  n(z)  =  1  +  6n(z)  where  6n  «  1. 

Typically,  6n«10  4  so  in  order  to  have  a  number  of  more  convenient  magnitude, 
the  " ref r activity"  N  has  been  defined  to  be  N  =  (n  -  1)  x  10®  =  fin  x  10^. 

Expanding  the  exponential  term  in  (8)  and  using  the  modified  refractiv- 
ity,  defined  as  M  -  N  +  10  z/a,  the  refractive  index  can  be  written  as 

n  =  1  +  10“6M  (9) 

Substituting  (9)  into  (7),  assuming  paraxial  rays  and  writing  the  angle  in 
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milliradians  (p  «  10  8)  yields  the  differential  forms 


“  •  t  i; dz '  =  10”3  f?  ** 


;io) 
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Integrating  the  first  expression  in  (10)  yields 


+  2 (K  -  M  ) 
o 


(ID 


where  p  is  the  angle  at  the  modified  refractivity  M  resulting  from  initial 

angle  R  at  M  Integrating  the  second  expression  in  (10)  with  the  assump- 
c  o. 

tion  that  dM/dz  is  constant  yields 


P  '  Pr 


X  = 


10  o 


(12) 


where  the  initial  condition  x  =  0,  p  «  6  is  assumed.  The  result  in  (11)  is 

o 

simply  a  stateme".  Snell's  law  in  the  paraxial  approximation  and  is  inde¬ 
pendent  of  the  vf4y  the  modified  refractivity  varies  with  height.  The  result 
in  (12),  no ever,  is  valid  only  if  the  modified  refractivity  varies  linearly 
with  height.  The  modified  refractivity  profile  for  the  clutter  model  is  thus 
assumed  to  be  linearly  segmented  with  height. 

Some  properties  of  ray  paths  which  are  important  to  the  question  of  sea 
clutter  undor  ducting  conditions  are  shown  in  Figure  2.  A  surface  based  duct 
resulting  from  an  elevated  refractive  layer  is  shown  diagramatically  in  Figure 
2  as  a  trilinear  M-profile.  Assume  that  a  ray  source  is  located  at  some 
height/  such  as  position  A  in  the  Figure.  From  (11)  a  ray  will  be  concave 
•upward  within  the  first  linear  segment  of  the  profile  and  concave  downward 
within  the  second.  Further/  the  magnitude  of  the  slope  of  any  ray  trajectory 
will  be  the  same  at  all  points  within  the  profile  with  the  same  M-value/ 
provided  of  course/  that  the  ray  reaches  the  particular  point  at  all.  For 
example,  a  ray  launched  horizontally  (labelled  [A]  in  Figure  2)  at  A  will  be 
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continually  turned  upward  until  reaching  the  first  change  in  profile  slope  and 
will  then  begin  turning  downward,  fis  the  ray  approaches  the  height  A' ,  where 
the  M-value  is  the  same  as  at  A,  the  ray  trajectory  will  again  approach  hori¬ 
zontal  and  continue  turning  downward,  tuvis  being  "reflected"  at  height  A' . 
The  ray  trajectory  again  changes  concavity  where  the  M- profile  slope  changes 
and  is  again  horizontal  at  height  A,  being  reflected  upward  at  this  height. 
The  ray  is  thus  "ducted"  or  "trapped"  between  A  and  A' . 

Consider  a  ray  launched  downward  at  A  with  such  an  angle  that  the  contin¬ 
ual  upward  turning  will  result  in  a  horizontal  trajectory  just  above  the 
surface.  Such  a  ray,  labeled  [B]  in  Figure  2,  will  be  trapped  between  the 
heights  of  B  and  B' ,  as  shown  in  Figure  2.  Any  ray  launched  below  this  ray 
will  reach  the  lower  surface  (which  is  necessary  before  the  ray  can  contribute 
to  sea  clutter)  and  be  reflected  according  to  the  usual  laws  of  reflection  at 
a  plane  interface.  There  is  a  launch  angle  for  which  a  ray  will  just  become 
horizontal  as  it  approached  C'  in  Figure  2.  Such  a  ray,  labeled  [C]  in  Figure 
2,  will  be  reflected  downward  at  C*  and  will  thus  be  trapped  between  the 
surface  and  C' .  Any  ray  launched  below  this  ray  (such  as  ray  [D]  in  Figure  2) 
will  reach  the  height  C'  at  a  non-zero  angle  and  continue  upward,  escaping 
through  the  top  of  the  duct.  There  is  thuf?  a  range  of  downward  launch  angles 
for  which  rays  will  be  both  trapped  in  the  duct  and  reflected  off  the  lower 
(sea)  surface.  By  symmetry  and  reciprocity  arguments,  there  is  also  a  family 
of  rays  with  upward  launch  angles  which  will  be  trapped  in  the  duct  as  well  as 
be  reflected  off  the  lower  surface.  It  is  these  two  families  of  rays  which 
contribute  to  sea  clutter  enhancement  under  conditions  when  an  elevated  re¬ 
fractive  layer  produces  a  surface  based  duct. 

In  the  ray-optical  propagation  model,  a  single  ray  is  envisioned  as 
departing  the  radar,  propagating  forward,  then  to  be  refracted  by  the  atmo- 


sphere  such  as  to  reach  the  surface  of  the  sea.  Some  of  the  radiated  energy 
is  returned  (on  a  reciprocal  path)  as  sea  clutter  while  t.he  bulk  of  the  energy 
is  reflected,  propagating  forward  to  be  again  refracted  by  the  atmosphere  to 
reach  the  surface  of  the  sea,  repeating  t’  .  procedure.  The  surface  back- 
ecatter  power  returned  to  the  radar  (sea  clutter)  can  be  conveniently  written 
in  terms  of  the  "radar  equation"  as 


G2  X2  2(m-1) 
pc  =  Pt — ~r  *e 

(4lt)3R4 

where  ?c  =  Clutter  power  returned 
Pt  -  Power  transmitted 
G  =  Antenna  power  gain 
K  ~  Radar  wavelength 

R  -  Radar- to- surface  distance  along  a  raypath 

r  =  Sea  reflection  coefficient 

m  =  Order  of  the  ray-hop 

a  =  Total  clutter  cross  section  at  range  R 
c 


(13) 


In  equation  (13),  the  sea  reflection  coefficient  is  taken  from  the  work  of 
Ament  (1953)  as  extended  by  Beard  (1961).  A  similar  expression  can  be  written 
for  a  target  by  replacing  the  term  r2^m  ^0  by  the  term  f40  (Kerr,  1951) 
where 

F  =  Pattern  -  propagation  factor  at  target  position 

o  =  Target  cross  section 
t 

The  definition  of  pattern-propagation  factors  is  such  that  the  ranges  at  which 
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target  detection  can  occur  (without  clutter)  are  given  by 

R  <  F  ®R,.  (14) 

is 

where  R^g  is  the  range  to  target  detection  in  free  space.  If  clutter  is 
present,  a  working  definition  of  target  detection  is  the  simultaneous  satis¬ 
faction  of  (14)  and  the  requirement  that  the  target  power  returned  exceed  the 
clutter  power  returned. 

R  <  F.Rfg  AND  PT  >  Pc  (15) 

where  PT  5  s  the  returned  target  power.  In  terms  of  the  various  cross 
sections,  (15)  is  equivalent  to: 

F  <  (  ac/Ct)  1/4r(m-1)/2  (16) 

This  model  is  based  on  a  ray-optical  propagation  model.  To  accomodate  a  non¬ 
ray  optics  model,  such  as  the  equivalent  single  mode  model  used  in  the  IREPS 
for  ducted  propagation  (Hattan,  1982),  note  tfc*fi  ,16)  can  be  used  provided 
only  that  a  way  to  determine  F  as  a  function  of  renne  is  available. 

The  area  extensive  cross-section  of  sea  return  is  a  function  of  many 
variables,  including  sea  state,  angle  of  incidence,  aspect  angle,  radar 
wavelength  and  polarization.  The  clutter  cross  section  model  adapted  for  use 
with  the  ray  path  popagation  model  was  developed  at  Georgia  Institute  of 
Technology  and  is  commonly  known  as  the  GIT  sea  clutter  model  (Trebits  et  al., 
1978).  The  model  is  semi-empirical  and  calculates  the  sea  clutter  cross 
section  normalized  per  unit  area  as  a  product  of  three  factors:  sea 


direction,  windspeed,  and  multipath.  Each  factor  is  a  function  of  the 
appropriate  independent  variables.  Thus,  or  =  0q/A,  where  cq  is  the 

normalized  cross  section  (see  Figure  3)  and  A  is  the  pulse  length  limited  area 
given  by 

A  =  ~  R  9 

where 

c  =*  speed  of  light 
x  =  pulse  length 
R  =  range 

0  =  horizontal  beam  width 

The  model  assumes  a  "fully  arisen"  sea  for  which  wave  height  and  w.^nd  speed 
are  coupled  under  the  assumption  of  equilibrium  conditions.  Both  the  wind/ -e a 
direction  dependence  and  the  wind  speed  dependence  of  the  cross  section  model 
are  emp  irically  derived.  The  multipath  portion  of  the  model  has  a  theoretical 
basis. 

The  sea  clutter  cross  section  per  unit  area  is  shown  in  Figure  3  for  two 
radar  frequencies  (C-Band  and  L-Band)  at  several  wind  speeds.  Both  upwind  and 
downwind  values  are  shown  for  a  15-knot  wind  speed.  Some  properties  of  the 
sea  clutter  cross  section  per  unit  are',  are  to  be  especially  noted.  For  high 
wind  speeds,  the  sea  clutter  cross  section  for  the  C-Band  increases  very 
rapidly  at  very  small  grazing  angles,  reaching  a  nearly  constant  value  between 
five  and  ten  milliradians.  This  property  is  also  exhibited  at  L-Band,  al¬ 
though  not  so  dramatically.  For  both  bands,  there  is  a  large  sensitivity  to 
variations  in  wind  speeds  for  moderate  speeds.  For  example,  there  is  as  much 
as  a  20  dB  change  in  cross  section  between  10-knot  and  15-knot  winds.  There 
is  a  much  decreased  sensitivity  to  changes  in  wind  speed  at  higher  wind 


speeds,  although  the  clutter  cross  sections  are  much  larger  for  higher  wind 
speeds.  Note  also  that  the  directional  dependence  of  clutter  cross  section  is 


insensitive  to  the  actual  windspeed  and  grazing  angle,  although  thers  is  some 
sensitivity  to  the  radar  band,  with  the  larger  directional  dependence  at 
higher  frequencies. 

The  computer  implementation  of  the  model  was  effected  in  the  following 
way.  Reference  modified  refractivities  are  computed  corresponding  to  the 
levels  in  Figure  2  at  C'  (denoted  by  Mu) ,  at  B  (denoted  by  Md)  and  at  A 
(denoted  by  Mj.) .  The  occurrence  of  a  surface  based  duct  is  determined  if  My 
<  Md  and  the  radar  is  in  the  duct  if  >  My.  If  these  are  satisfied,  the 
launch  angle  ( Pd)  for  the  surface  grazing  ray  ( [B]  in  Figure  2)  and  the  launch 
angle  (  p^)  for  the  ray  just  grazing  at  C' ( [C]  in  Figure  2)  are  computed.  The 
rays  that  are  trapped  and  also  contribute  to  sea  clutter  are  found  within  the 
limits  I  Pu  "  Pdl  both  above  and  below  the  horizontal.  This  fan  of  rays  is 
divided  into  30  equally  spaced  rays  ( p^ ) ,  and  for  each  ray  the  following  are 
computed: 

P  .  =  grazing  angle  at  sea  surface  for  j-ray. 

SD 


=  surface  distance  for  the  upward  launched  j-ray  to  be  "reflected" 
and  return  to  the  sea  surface. 


surface  distance  for  the  downward  launched  j-ray  to  re*  ch  the  sea 
surface. 


Ap  =  skip  distance  for  the  j-ray  =  0  +  0 

3  *uj  Kdj 


0  ,  =  sea  clutter  cross  section  for  j-ray  (normalized  per  unit  area) 
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r.  =  sea  surface  reflection  coefficient  for  j-ray 
3 

The  p  are  computed  using  (11)  and  the  p's  are  computed  by  iterative  applies- 
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tion  of  (12)  for  each  ray. 

The  propagation  path  is  divided  into  100  range  bins  of  2  ksa  each.  A 
recursive  procedure  is  then  employed  to  determine  the  contribution  to  the 
total  clutter  power  due  to  each  successive  hop  for  both  the  upward  launched 
and  downward  launched  ray  sets.  Values  within  any  given  range  bin  are  deter¬ 
mined  by  linear  interpolation  between  two  j-rays  if  a  single  ray  does  not 
"land"  within  the  bin.  The  100  elements  of  the  range  bin  array  are  finally 
modified  to  include  a  preset  but  arbitrary  target  cross  section  and  appropri¬ 
ate  roots  are  taken  to  provide  the  right  hand  side  of  (16). 

The  operational  implications  of  the  combined  ray  path  propagation  model 
and  normalized  clutter  cross  section  model  can  be  qualitatively  inferred  from 
the  considerations  of  the  individual  aspects  already  discussed.  It  can  be 
expected,  for  example,  that  at  high  wind  speeds,  even  moderately  strong 
surface  based  ducts  will  result  in  enhanced  clutter  effects.  The  clutter  ring 
"splitting,"  mentioned  previously  (Snyder,  1979)  in  regards  to  the  USS  ENTER* 
PRISE  demonstrations  of  the  IREPS,  can  be  inferred  from  the  existence  of  two 
sets  of  influential  rays  (the  initially  upgoing  and  downgoing  launched 
sets) .  The  existence  of  clutter  rings  in  only  a  limited  angular  sector  might 
be  inferred  by  the  directional  dependence  of  the  clutter  :ross  section,  but 
horizontal  inhomogeneity  in  the  environment  would  seem  a  more  likely  cause. 

The  clutter  effects  model  was  coded  and  implemented  in*o  a  special  FOR¬ 
TRAN  version  of  the  IREPS.  Typical  results  are  presented  in  the  "coverage 
diagrams"  of  Figures  4-6.  A  hypothetical  C-Band  radar  with  100-nmi  free  space 
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detection  range  against  a  reference  target  of  two  square  metres  cross  section 
is  assumed  located  at  140  feet  height.  As  a  reference,  the  coverage  diagram 
for  a  nonducted  (standard)  atmosphere  with  a  10-knot  wind  is  shown  in 
Figure  4.  Recall  '-.:';iat  the  darkened  area  of  the  diagram  represents  the 
range/height  conditions  where  the  radar  would  detect  the  target. 

A  coverage  diagram  for  the  same  wind  conditions,  but  for  a  surface  based 
duct,  is  shown  in  Figure  5.  The  elevated  refractive  layer  causing  this  duct 
is  located  at  about  300  metres  and  is  40  metres  thick.  The  M-unit  deficit 
between  the  surface  and  the  layer  minimum  is  40  M-units.  Notice  the  enhanced 
coverage  for  low  altitudes  to  longer  ranges.  For  this  case,  any  clutter 
effects  are  just  barely  observed.  A  coverage  diagram  for  the  same 
refractivity  profile  but  for  a  15— knot  wind  is  shown  in  Figure  6.  For  this 
case,  tie  predicted  effects  of  clutter  are  quite  dramatic.  Within  the  darken¬ 
ed  areas  of  the  diagram,  target  detection  would  be  expected.  However,  within 
the  areas  no  longer  darkened  compared  to  Figure  5,  the  target  would  be  unde¬ 
tectable  compared  to  clutter.  The  presence  of  clutter  in  this  case  would 
appear  on  a  PPI  type  display  as  bright  rings,  masking  the  target. 


EFFECTS  OF  REFRACTIVE  LAYER  BOUNDARY  PERTURBATIONS 


An  important  property  of  the  clutter  model  is  the  ability  to  predict 
characteristics  of  clutter  rings.  A  crucial  factor  in  this  ability  is  the 
properties  of  the  vertical  modified  refractivity  profile.  Variations  in  the 
profile  boundaries  can  be  expected  to  produce  variations  in  the  clutter  ring 
characteristics.  For  this  reason,  effects  of  both  bulk  and  localized 
variations  in  profile  boundaries  have  been  examined. 

Typical  sensitivities  of  computed  clutter  ring  locations  to  bulk 

variations  in  boundary  heights  for  elevated  layers  are  shown  in  Figure  7.  The 

figure  shows  the  location  of  the  first  and  second  clutter  rings  for  a 

trilinear  profile  as  a  function  of  layer  thickness  and  parametric  in  layer 

height.  In  this  example,  variation  in  layer  thickness  results  from  changes  in 

the  slope  of  the  modified  refractivity  within  the  elevated  layer  with  no 

change  in  the  height  of  the  bottom  of  the  layer,  the  slope  of  the  lowest 

region  of  the  profile  (called  AM/ Ah  in  Figure  7),  or  the  net  M-deficit 

o 

(called  Am  in  Figure  7).  The  geometrical  change  in  the  profile  is  shown  in 
the  figure  insert.  The  change  in  clutter  ring  location  per  metre  change  in 
layer  thickness  is  approximately  0.2  nmi  for  the  first  ring  and  0.3  nmi  for 
the  second.  These  results  depend  only  very  slightly  on  the  layer  height  or 
thickness. 

An  alternate  variation  in  layer  thickness,  resulting  from  raising  the 
height  of  the  layer  while  maintaining  a  fixed  gradient  both  within  and  below 
the  layer  as  well  as  a  fixed  M-deficit  was  also  considered.  For  this  situa¬ 
tion,  the  change  in  clutter  ring  location  per  metre  change  in  layer  height  is 
approximately  0.1  nmi  for  the  first  ring  and  0.2  nmi  for  the  second.  Just  as 
for  the  previous  case,  these  results  are  insensitive  to  the  specific  layer 
height  or  thickness. 
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From  these  and  other  similar  results/  it  is  concluded  that  bulk  changes 
in  refractive  layer  height  or  thickness  on  the  order  of  a  few  metres  will  only 
slightly  change  the  apparent  clutter  ring  location  on  a  typical  PPI  display. 
Bulk  changes  of  such  an  order  might  be  expected  over  short  periods,  perhaps  on 
the  order  of  a  few  tens  of  minutes.  Over  periods  of  a  few  hours  or  longer, 
bulk  changes  on  the  order  of  tens  of  metres  in  refractive  layer  height  or 
thickness  may  be  expected.  For  such  longer  periods,  changes  of  a  few  miles  in 
clutter  ring  location  might  be  expected. 

The  actual  occurrence  of  clutter  rings  is  a  more  sensitive  function  of 
refractive  layer  configuration  than  is  the  variation  in  location.  Changes  in 
layer  height  or  thickness  are  usually  accompanied  by  changes  in  M-deficit. 
This  results  in  changes  in  the  grazing  angle  and  consequently  in  the  clutter 
return  signal.  Bulk  changes  in  height  or  thickness  of  elevated  refractive 
layers  are  thus  expected  to  lead  to  clutter  rings  which  fade  "in-and-out"  on  a 
typical  PPI  display  while  remaining  nearly  fixed  in  location. 

The  effects  of  bulk  changes  in  the  boundaries  of  the  refractive  layer 
formed  by  a  surface  evaporation  layer  have  also  been  considered.  The 
evaporation  layer  was  modelled  simply  with  a  single  negative  M  gradient  at  the 
surface,  changing  to  a  positive  gradient  at  height  hQ  as  shown  in  the  insert 
to  Figure  3.  Although  ray  optical  methods  are  generally  inadequate  for 
complete  quantitative  analysis  of  propagation  in  the  evaporative  duct,  the 
method  is  believed  adequate  to  describe  qualitative  and  quantitative  trends 
resulting  from  modifications  in  evaporative  duct  parameters.  Two  types  of 
changes  were  considered,  one  with  constant  layer  strength  (height)  and  vari¬ 
able  M-deficit  and  the  other  with  fixed  M-deficit  and  variable  layer 
strength.  Typical  results  for  variable  M-deficit  are  shown  in  Figure  8.  In 
the  figure,  an  evaporation  layer  with  minimum  H-value  at  ten  metres  and  vari- 
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able  K-deficit  is  superimposed  on  a  modified  refractivity  profile  with  a 
50-matre-thick  refractive  layer  at  a  height  of  300  metres.  The  original 
profile  (without  evaporation  layer)  produces  a  ground  based  duct  with  a  30  M- 
unit.  deficit.  A  diagramatic  picture  of  the  modified  refractivity  profile  is 
included  in  Figure  8. 

Figure  8  shows  that  the  location  of  clutter  rings  is  very  insensitive  to 
the  change  in  M-deficit.  In  fact,  for  the  range  of  M-def icits  considered,  the 
change  in  clutter  ring  location  is  less  than  0.1  nmi.  However,  the  clutter 
return  cross  section  is  seen  to  be  a  sensitive  function  of  the  M-deficit, 
increasing  in  excess  of  3  dfi  for  the  range  in  M-daficit  shown.  Thus,  for  this 
type  of  bulk  profile  modification,  clutter  rings  would  exhibit  exceptional 
spatial  stability  with  variable  intensity. 

The  alternate  variation  in  evaporation  layer,  characterized  by  a  fixed 
value  to  M-deficit  and  variable  layer  height  (hQ  of  Figure  8),  showed  a  signi¬ 
ficant  sensitivity  of  computed  ring  location  to  layer  height  (or  layer 
strength) .  Because  the  M-deficit  remained  constant,  and  therefore  also  the 
grazing  angle,  the  clutter  model  predicts  this  type  of  profile  modification 
produces  no  change  in  the  clutter  return  signal. 

The  variations  in  clutter  ring  characteristics  dm  to  modelled  variations 
of  an  evaporation  layer  thus  differ  both  qualitatively  and  quantitatively  from 
the  results  for  variations  of  an  elevated  layer.  In  particular,  clutter  ring 
location  is  reasonably  stable  for  bulk  variations  in  elevated  layers,  but 

highly  variable  for  bulk  variations  in  evaporation  laye  .'•>  Furthermore,  much 
larger  changes  in  ring  intensity  would  be  expected  for  evaporation  layer 

variations . 

Localized  variations  in  layer  boundaries  were  modelled  in  terms  of  spa¬ 
tially  periodic  fluctuations  of  layer  boundaries.  Two  forms  of  perturbations 
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were  considered/  one  a  sinusoidal  variation  in  evaporation  layer  strength 
(height)  and  the  other  an  exponentially  damped  sinusoidal  variation  in  refrac- 
tivity  at  the  boundary  of  an  elevated  refractive  layer.  Numerical  integration 
of  the  three-dimensional  ray  equation  was  used.  The  eikonal  equation  in 
rectangular  component  form  (see  equation  5)  was  employed  with  the  assumption 
that  9n/9y  =  0  so  that  n  had  only  a  two-dimensional  variation  ( z-vertically ; 
x-in  the  horizontal  propagation  direction) .  Also,  an  arbitrary  initial  ray 
direction  and  position  were  permitted.  The  coupled  equations  (5)  were  solved 
with  a  fourth-order  Runge-Kutta  routine  using  a  comparison  with  second-order 
to  control  step  size  adjustment. 

The  use  of  the  raytrace  formulation  of  (5)  requires  an  analytic  expres¬ 
sion  for  the  refractive  index.  In  order  to  approximate  the  trilinear  profiles 
used  previously,  the  procedure  of  Booker  (1977)  was  adopted.  The  steps  in 
this  procedure  are  depicted  in  Figure  9.  A  prototype  modified  refractivity 
profile  composed  of  linear  segments  with  slopes  A^  is  assum  i.  The  transition 
from  one  slope  to  another  is  modelled  by  the  analytic  expression  (A) .  By  an 
appropriate  selection  of  the  "scale  height"  (H),  the  slope  rapidly  approaches 
A^  for  z  <  z2  and  A2  for  z  >  z2>  A  generalization  to  a  multislope  profile  is 
given  by  expression  (B).  The  final  analytic  form  of  the  modified  refractivity 
is  found  by  integrating  expression  (B)  and  is  given  by  expression  (C).  A 
sample  of  the  use  of  the  algorithm  to  generate  profiles  is  shown  in  Figure 
10.  In  the  figure,  the  prototype  profile  is  trilinear  and  the  results  for  two 
scale  heights  (H  =  20  km  and  5  km)  are  included.  A  scale  height  of  5  km  was 
used  for  the  modified  refractivity  profiles. 

Numerous  raytraces  were  computed  using  a  variety  of  modified  refractivity 
profiles.  Sample  raytraces  showing  the  effects  of  the  sinusoidal  pertur¬ 
bations  are  given  in  Figure  11.  The  effects  of  a  periodic  elevated  layer  are 
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shown  in  Figure  11a.  The  unperturbed  refractivity  profile  is  the  same  as  in 
Figure  10.  The  sinusoidally  perturbed  modified  refractivity  is  given  by 

M(x,z)  =  M  (z)  [1  +  A  exp((z  -  z)2/\2)  sin(2lcx/X  )] 
o  o  z  x 

where  MQ(z)  is  the  unperturbed  refractivity,  A  =  10,  zQ  =  300  m,  X^  =  5  m, 
and  X^  =  10  km.  The  significant  point  of  Figure  11a  is  the  fact  that  the 
range  to  reflection  or  turning  points  of  the  lowest  order  hops  is  only 
slightly  affected  by  the  perturbation.  What  is  significantly  affected  is 
whether  or  not  a  ray  reaches  the  surface.  These  characteristics,  which  were 
found  to  be  typical  in  general  for  the  elevated  layer  refractivity  perturba¬ 
tions  considered,  would  result  in  clutter  rings  "fading  in-and-out"  with  only 
slight  fluctuations  in  location. 

The  effects  of  a  periodic  perturbation  to  the  evaporation  layer  are 
shown  in  Figure  11b.  The  unperturbed  profile  is  the  same  as  for  Figure  11a 
but  with  a  model  evaporation  layer  such  as  shown  in  Figure  8  appended.  in 
this  case,  the  layer  strength  (height,  hQ)  is  20  m  and  the  additional  M- 
deficit  is  20  M-units.  The  sinusoidally  perturbed  layer  height  is  given  by 

h(x)  =»  hQ  +  A  sin(2nx/Xx) 

where  hQ  is  the  unperturbed  height  (20  m) ,  A  =  2  m,  X^  =  2  knu  Most  of  the 

"'ays  of  Figure  11b  are  only  slightly  modified  by  the  layer  perturbation. 
However,  a  significant  effect  is  demonstrated  on  one  ray,  where  it  is  seen 

that  the  perturbation  causes  trapping  within  the  evaporation  layer.  This 
trapping,  which  would  result  in  an  increased  clutter  return  from  an  extended 
area,  is  very  sensitive  to  the  properties  of  the  evaporation  layer. 
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Just  as  for  the  bulk  profile  modifications,  the  variations  in  clutter 
ring  characteristics  differ  for  evaporation  layer  fluctuations  as  compared  to 
elevated  layer  fluctuations.  Clutter  ring  locations  are  very  variable  for 
evaporation  layer  fluctuations  but  quite  stable  for  elevated  layer  fluctua¬ 
tions.  These  characteristics  are  the  same  as  for  bulk  variatians  in  refrac¬ 
tive  layers.  Thus,  clutter  signal  characteristics  on  a  typical  PPI  display 
might  indicate  the  type  of  refractive  layer  fluctuations,  whether  evaporation 
layer  or  elevated  layer.  The  characteristics  could  provide  only  limited 
information  about  the  extent  of  the  fluctuations,  whether  localized  or  wide- 


MODEL  VALIDATION 


The  sea  clutter  model  presented  here  employs  a  union  of  classical  one¬ 
dimensional  ray-optics  propagation  with  a  semiempirical  normalized  sea  clutter 
cross  section  model  not  originally  intended  for  application  in  this  manner. 
The  model  predictions  are  in  agreement  with  what  limited  data  could  be  found 
for  sea  clutter  return  under  ducting  conditions.  These  data,  in  general, 
concern  the  geometrical  aspects  and  not  the  amplitude  aspects  of  the  clutter 
signals.  Therefore,  it  is  recognized  that  some  sort  of  quantitative 
"validation"  effort  should  be  conducted  before  implemention  into  the  IREPS  for 
shipboard  use. 

Because  of  manpower  and  money  limitation;?,  an  extensive  validation  effort 
has  not  been  possible.  However,  a  NOSC  tenant  activity,  the  Integrated  Combat 
Systems  Test  Facility  (ICSTF),  agreed  to  permit  use  of  one  of  the  test  radars 
on  a  limited,  "not-to- interfere"  basis  for  the  acquisition  of  model  validation 
data:  The  radar  employed  was  an  operational  fleet  type  C-band  radar  (an  SPS- 
10)  permanently  shore  base  mounted  overlooking  the  SQCAL  OPAREA  (Southern 
California  Operating  Area) . 

Modification  of  the  SPS-10  necessary  to  permit  sampling  of  return  signals 
is  diagrammatic ally  depicted  in  Figure  12.  Data  are  recorded  on  cassette  tape 
using  a  microprocessor  controlled  Data  Acquistion  System  (DACQ) .  The  micro¬ 
processor  permits  sampling  with  a  two-minute  averaging  of  as  many  as  100  range 
bins  (spaced  from  0.25  nmi  to  1.22  nmi)  at  each  of  18  or  fewer  directions  with 
a  variable  angular  spacing  of  from  approximately  three  degrees  to  eleven 
degrees.  The  fundamental  timing  for  the  DACQ  is  derived  from  the  SPS-10 

trigger.  The  target  return  signal  is  sampled  at  the  IF-Amplifier,  prior  to 
detector  circuits,  and  passed  through  a  log- de ten cor  before  averaging  and 


recording.  Provisions  for  display  of  data  sample  point  locations  on  a  PPI 


display  unit/  along  with  conventional  PPI  target  display,  are  included. 

Some  preliminary  data  were  obtained  in  order  to  check  out  the  operation 
of  the  DACQ,  particularly  the  timing  operations.  One  set  of  data  was  obtained 
in  January  1980  during  a  period  of  moderately  high  winds  \\V.  Xttc  average  with 
gusts  to  20  kts) .  Sea  clutter  return  signals  for  this  case  are  shown  in 
Figure  13.  Shown  in  the  Figure  are  data  for  two  different  two-m.1'  -'he 
averages,  the  upper  corresponding  to  18  radials,  spaced  6.3°  apart  for  &  total 
of  107°  with  range  cells  every  0.97  nmi,  and  the  lower  corresponding  to  18 
radials  spaced  3.5°  apart  for  a  total  of  60°  with  range  cells  every  0-32 
nmi.  The  two  sets  of  curves  have  nearly  identical  range  dependence,  indica¬ 
ting  proper  operation  of  the  DACQ  timing.  The  large  signal  centered  around 
18.5  nmi  in  the  upper  figure  is  from  the  Coronado  Islands,  south  of  San  Diego. 

A  second  set  of  data  was  obtained  in  April  during  a  period  of  ducted 
propagation  but  very  low  wind  speeds.  As  before,  two  different  range  bin  and 
angular  separations  were  employed.  A  plot  of  the  radar  return  signal  is  shown 
in  Figure  14.  The  keypoint  in  the  figure  is  the  strong  signal  from  about  55 
nmi  to  65  nmi,  which  is  the  correct  range  for  the  Santa  Catalina  and  San 
Clemente  Islands.  Figure  15  shows  an  additional  verification  of  this.  The 
figure  is  a  map  of  a  portion  of  the  SOCAL  OPAREA.  The  straight  lines  from  San 
Diego  are  the  data  bearings,  according  to  the  DACQ  timing,  and  the  darkened 
areas  on  the  lines  represent  the  radar  return  signals.  The  positions  of  San 
Clemente  and  Santa  Catalina  Islands  are  clearly  indicated  in  the  signal 
return.  The  return  on  the  central  path  corresponds  to  the  location  of  the 
very  small  Santa  Barbara  Island,  which  is  not  shown  on  the  map.  The  "Xs" 
indicate  positions  of  ships  reporting  sea  state  information  at  any  time  during 
the  three  hours  preceding  or  following  the  data  recording.  The  first  numeral 
represents  the  sea  wave  height  and  ^he  second  represents  the  swell.  Note  that 
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the  waves  were  not  sufficient  to  produce  clutter  signals. 

Based  on  these  and  other  similar  data,  it  was  concluded  that  the  DACQ 
unit  operated  properly.  The  model  validation  requires  not  just  the  occurrence 
of  ducted  propagation,  which  happens  frequently  in  the  SOCAL  area,  but  also 
the  occurrence  of  a  sea  state  sufficient  to  provide  significant  clutter 
return.  The  occurrence  of  both  together  is  very  rare a  The  procedure  followed 
was  to  anticipate  the  occurrence  of  ducted  propagation,  such  as  through 
analysis  of  local  weather  patterns,  and  to  operate  the  DACQ  whenever  such 
propagation  conditions  were  available,  provided  the  use  of  the  SPS-10  would 
not  interfere  with  dedicated  ICSTF  activities.  The  limitation  to  normal 
working  hours,  which  automatically  eliminated  weekend  recording  periods,  the 
weekly  service  schedule,  which  eliminated  one  day  a  week,  and  the  "normal" 
radar  usage,  which  averaged  about  two  days  a  week,  left  little  time  for  a 
"not-to- interfere"  usage.  Nevertheless,  several  hours  of  return  signals  wore 
recorded  during  ducted  propagation  conditions.  Unfortunately,  only  one  set  of 
data,  obtained  in  September  1980,  had  a  sea  state  sufficient  to  result  in  sea 
clutter  return. 

The  propagation  conditions  occurring  during  ♦-his  time  are  shown  in  Figure 
16.  Figure  16a  shows  refractivity  condi tons  obtained  in  the  San  Diego  area  on 
1  'ptember  30,  1930  at  1200  2  and  2155  Z.  The  refractivity  layer  formed  at  an 
altitude  of  approximately  300  metres  in  the  near  shore  region  and,  as  can  be 
seen,  was  very  stable  for  at  least  a  twelve-hour  period  with  an  M-deficit  of 
approximately  30  M-units.  Figure  16b  shows  refractivity  profiles  obtained  at 
Point  Mugu,  north  of  San  Diego,  for  the  same  time  period.  In  this  case  the 

refractivity  laytr,  also  formed  around  300  metres,  was  more  variable,  but 
nevertheless  persisted  over  the  full  twelve-hour  time  period.  The  M-deficit 
for  Point  Munu  varied  from  a  value  in  exces  of  40  M-units  to  as  small  as  15 


M-'inits  for  this  time  period.  The  wind  speed  at  surface  level  reported  by 
ships  in  the  area  varied  from  slight  to  moderate.  Visual  observations  at  NOSC 
indicated  moderate  white  cap  formation. 

A  normalized  sea  clutter  power  versus  range  was  calculated  using  the 
models  previously  discussed.  An  example  of  thene  results  is  shown  in  Figure 
17  for  the  refractivity  profiles  obtained  at  NOSC  at  2155  Z.  The  key  point  to 
notice  is  the  large  increase  in  clutter  power  at  a  range  of  approximately 
40nmi  and  the  second  increase  at  a  range  ol‘  approximately  90  nmi.  The  second 
increase,  however,  only  reaches  a  level  approximately  15  dB  below  the  40~nmi 
results.  The  clutter  power  results  presented  in  Figure  17  are  quite  typical 
of  the  results  obtained  for  the  profiles  persisting  throughout  the  time  period 
presented  in  Figure  16. 

Radar  backscatter  signals  averaged  over  a  two-minute  period  were  recorded 
at  approximately  30-minute  intervals  during  the  same  time  period  of  Figure 
16.  Results  obtained  for  a  short  period  between  2000  Z  and  2300  Z  are  shown 
in  Figure  18.  The  presentation  in  Figure  18  shows  essentially  a  tracing  of 
what  would  be  displayed  on  a  PPI  radar  scope,  superimposed  on  a  map  of  the 
geographic  area.  Note  the  darkened  outline  around  San  Clemente  Island.  This 
represents  the  large  signals  returned  from  the  island.  The  signals  returned 
from  very  large  distances,  such  as  from  Santa  Barbara  Island,  which  is 
approximately  100  nmi  from  the  radar,  indicates  the  generally  extensive 
persistence  of  the  surface  based  ducting  region  which  occurred  during  this 
time  period.  Because  only  limited  angular  sectors  could  be  examined  at  any 
given  time  with  the  data  acquisition  system,  not  all  angular  sectors  are 
included  in  this  time  period.  For  those  sectors  which  are  included  there  is 
generally  a  significant  signal  return  from  the  30-  to  50-nmi  range,  which  is 
consistent  with  the  model  predictions.  The  relative  signal  strengths  from  the 
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various  sectors  differed  considerably,  by  as  much  as  15  to  20  <3B.  Because 
each  sector  was  sampled  au  different  times,  these  variations  are  possibly 
caused  by  spatial  and  temporal  variations  in  refractivity  and  wind  speed. 

A  major  requirement  in  the  validation  of  the  clutter  model  is  the  avail¬ 
ability  of  the  clutter  return  absolute  amplitude.  The  measurement  of  appro¬ 
priate  voltages  along  the  radar  input  circuits  while  using  a  calibrated  input 
at  the  antenna  was  not  possible  due  to  radar  use  limitations.  The  calibration 
procedure  used  employed  a  small  calibrated  target  (of  known  cross  section) 
being  towed  radially  away  from  the  radar  site.  The  target  was  mounted  on  a 
2.5-metre  mast,  made  of  plastic  sewer  pipe,  and  the  mast  was  mounted  on  a  » 
small  rubber  raft  using  all  nonconducting  materials.  The  raft  and  mast  ar¬ 
rangement  was  then  towed  at  slow  speed  (5  kts)  using  an  800-foot  tow  line.  In 
the  intervening  time  between  the  acquisition  of  the  September  30,  1980  data 
and  the  calibration  data,  numerous  modifications  were  made  in  the  SPS-10  radar 
by  ICSTF  technicians.  Included  were  changes  in  the  antenna  feeds  as  well  as 
numerous  electronic  components  within  the  radar  proper.  Because  of  this, 
absolute  calibration  for  the  September  1980  data  i3  unavailable. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A  radar  sea  clutter  model  based  on  a  ray-optics  propagation  formulation 
combined  with  a  semi-empirical  normalized  sea  clutter  radar  cross-section  has 
been  presented.  A  preliminary  form  of  the  model  was  coded  in  FORTRAN  for 
incorporation  into  the  IREPS.  The  model  was  exercised  for  typical  propagation 
conditions  and  round  to  be  in  agreement  with  what  limited  observation  data  are 
available.  A  single  frequency,  single  location  validation  effort  was  conduct¬ 
ed.  Model  and  observed  results  were  in  agreement,  although  absolute  amplitude 
data  were  not  obtained. 

A  multifrequency  and  multilocation  validation  effort  should  be  conducted 
using  dedicated  equipments.  The  locations  should  be  selected  to  optimize  the 
occurrence  of  ducted  propagation  and  significant  sea  clutter  return. 
Incorporation  of  clutter  models  into  IREPS  should  be  done  only  after  a 


successful  validation  effort 
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Figure  1.  Transformation  from  cylindrical  coordinate 
frame  to  Cartesian  coordinate  frame. 
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Figure  2.  Sample  ray  paths  with  a  surface  based  duct. 


HEIGHT  IN  FEET  HEIGHT  IN  FEET  HEIGHT  IN  FEET 


Figure  4.  Coverage  diiigram  for  hypothetical  C-band  radar  - 
standard  (non  ducted)  atmosphere. 


Figure  5.  Coverage  diagram  for  hypothetical  C-band  radar 
surface  based  duct.  Duct  height  -  300  metres;  refractive 
layer  thickness  -  30  metres. 
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Figure  6.  Coverage  diagram  for  hypothetical  C-band  radar  - 
surface  based  duct.  Duct  height  *  300  metres;  refractive 
layer  thickness  -  30  metres. 
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Figure  10.  Sample  of  analytic  profile  generation. 
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Figure  1  lb.  Sample  effect  of  periodic  evaporation  layer. 
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Figure  12.  Modifications  to  SPS-10  radar. 
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Figure  13.  Sample  signals  with  nonducted  propagation. 
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Figure  14.  Sample  signals  with  ducted  propagation. 


Figure  1 5.  Map  of  SOCAL  OP-AREA  depicting  radar  return. 
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Figure  ’6.  Modified  refi  activity  profiles 

A)  San  Diego,  California 

B)  Point  Mugu,  California. 
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Figure  17.  Normalized  clutter  power  (computed). 
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